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(57) Abstract: A method for conversion of input audio frequency data, at an input sample frequency, to output audio fr«]uency data, 
at an output sample frequency. The input data is subjected to expansion to produce expanded data at an output sample frequency. The 
expanded data is interpolated to produce output data. In one embodiment of the invention the interpolation is effected by a process 
that also filters the ou^ut data. In another embodiment, the input data is sampled by an integer factor to produce expanded data, 
the expanded data is then interpolated to produce the ou^ut data. Also disclosed is a method of transition of a signal output, at one 
frequency, to a signal output al another fiequency. The signal output at said one frequency is faded out over a period, and the signal 
output at said other frequency is faded in over that period. Both signal outputs are combined to produce the signal output over said 
period. Apparatus for effecting the methods is also disclosed. 
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UNIVERSAL SAMPLING RATE CONVERTER FOR DIGITAL 
AUDIO FREQUENCIES 



5 FIELD OF THE INVENTION 

This invention relates to methods and apparatus for conversion of input audio 
ftequencies to output audio fiequencies, by expanding input data and interpolating the e3q>anded 
data to form the output data at a new frequency. 

10 

BACKGROUND OF THE INVENTION 

Distal Audio is based on many different means of commnnication, as described in the 
following reference. The diflferent digital media generally have conflicting sampling fiequendes, 
15 For example, digital transmission of broadcasting program at 32 kHz, compact discs at 44.1 
kHz, digital video discs at 48 kHz and speech recording at 6 kHz to 8 kHz, as described in 
"High Quality Digital Audio in the Entertainment Industry", IEEE ASSP Magazine, 1985 pages 
2-25. Thus, digital audio requires a sampling frequency conversion technique to handle simple 
as well as non-trivial ratios efficiently. 

20 

Conversion by going from digital to analogue (through a DAC and a low-pass filter) and 
then re-sampling the smoothed signal at the ouQ)ut rate is simple, but costly and limited by the 
imperfections (non-linearity, phase response, noise) of the analogue filter as desmbed in **HigJi 
Quality Analogue Filters for Digital Audio", 67* AES Convention, November 1980. 

25 

Conversion in simple integer or rational ratios f/f^ by single or multi-stage FIR filter 
design, as described in Rabiner and Croichiere, Multi-rate Digital Signal Processing, Prentice 
Hall Publication, 1983. However, it is not particularly suited for many arbitrary ratios, as it leads 
to far too many filter configurations. An individual filter configuration is suited maximally to 
30 a subset of these ratios only. 
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SUMMARY OF THE INVENTION 

According to the processes of the invention there is provided a method for conversion 
5 of input audio frequency data, at an input sample frequency, to output audio frequency data, at 
an output sample frequency. The method including the steps of subjecting the input data to 
expansion to produce expanded data at an output data sample frequency; and interpolating the 
expanded data to produce the output data, wherein the interpolation process also filters the 
output data. 

10 

According to the processes of the invention there is also provided a method for 
conversion of digital input audio frequency data, at an input sample frequency, to digital output 
audio frequency data, at an output sample frequency. The method including upsampling the 
input data by an integer factor, so as to increase the sampling rate of the ii3put data to produce 
1 5 expanded data; and interpolating the expanded data to produce the oxitput data. 

The invention includes the method of transition of a signal output at one frequency to 
signal output at another frequency. In this method the signal output at said one frequency is 
faded out over a period, and the signal output at said another frequency is faded in over that 
20 period, and both are combined to produce the signal output over said paiod. 

In accordance with the invention an apparatus for implementing transition of a signal 
output thereof from a condition at which the signal output is derived from a first input signal, 
of a first frequency, applied on a first input to the apparatus, to a condition at vs^ch the signal 
25 output is derived fix>m a second input signal, of a second frequenqr, applied to a second input 
of apparatus, the apparatus having means for fading out th& first input signal over a period and 
fading in the second input signal over said period, and means for combining these over said 
period and applying the so combined input signals to said signal output. 

30 According to the invention ttiere is provided a digital sequential frequency converter for 

conversion of input audio frequency data, at an mput sample frequency, to output audio 
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frequency data, at an output sample frequency. Included in the apparatus is the means for 
subjecting the input data to expansion to produce expanded data at an output data sample 
frequency; and the means for interpolating the expanded data to produce the output data, by an 
interpolation process which also filters the ou^ut data. 

5 

According to the invention there is also provided a digital sequential frequency converter 
for conversion of digital input audio frequency data, at an input sample frequency, to digital 
output audio frequency data at an output sample frequency. Included in the apparatus is the 
means for iq)sampling the input data by an uatego: fector, so as to increase samplii^ rate of the 
10 input data to produce expanded data; and the means for intapolating the expanded data to 
produce output data. 

A single simple stractiire is often desired in Audio applications, for conversion between 
commonly occurring frequencies. The advantage of using a single structure is that for 
15 conversion between different frequency combinations, the same block code and same 
coefficients can be used. This reduces the program code size. A single simple structure also 
means it can be implemented efficiCTtly as a hardware block, without excessive chip area. 

A single and simple embodiment of the invention may be implemented on a two-stage 
20 design for sampling rate conversion for a group of audio frequencies. Embo dime n t s of the 
invention may be implemented in this way. 

Grenerally, the samplii^ and reconstruction of data may be effected in accordance 
with a suitable reconstruction foimxila. Thus, Let x[n] be sample of band-limited signal x(t). 
25 with sampling frequency of Fs. Then samples x'[m] corresponding to sampling frequency F' 
can be computed using the reconstruction formula: 
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x^[m]=Z :c[n + fc]<()jt(Ar) where 



and 0)^ = min (tcF^, ^F'). Thra n = INT [mT*/TJ and At = mT' - nTs. 

A two-stage embodiment may also be employed In the first stage an npscaling by a 
5 factor of 1 6 may be effected, by means of zero insertion (fifteen zeros for every one sample) 
followed by fiiltering with a fiOiter cut-off firequency of tz/16. The zero insertion and the filtering 
may be performed in practice through polyphase method that considerably decreases the 
computation complexity. 

10 The second stage may be a linear or sine interpolator. In simple cases, linear 

intCTpolation may be used. In that linear interpolation is used a fijrther reduction in Hxe 
computation can be achieved by ignoring those i)olyphase outputs which do not contribute to 
the actual output. 

15 InUxe case where fixed buffer constraint exists and slight variation occurs in input and 

ou^ut sampling clocks, method of cross fading may be used to smoothen the output. 



DETAILED DESCRIPTION OF THE DRAWINGS 

20 

The invention is further described by way of example only with reference to the 
accompanying drawings, in which: 



Figure 1 is a block diagram of a digital fiequency converter of a general kind described 
25 as prior art; 
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Figure 
with the invention; 



2 is a block diagram of a digital frequency converter construction in accordance 



Figure 3 is a block diagram of a digital frequency converter construction in accordance 
5 with the invention; 

Figure 4 is ablock diagram of a digital fi^quency converter constmctio^ 
with the invention; 

10 . FigureSisaflov^diagramdepictingsincandlinearinterpolationte^^^ 
be efEected in the figure 3; 

Figure 6 is a diagram depicting processing steps occurring in the use of the digital 
fiequencq' converter of figure 3; and 

^ ^ Figure 7 is a diagram illustrating a cross-feding technique, in accordance with the 

invention. 



20 



DETAILED DESCRimON OF THE PREFERRED EMBODIMENT 



Figure 1 shows an audio ftequency converter 1 0 according to the prior art Hiis employs 
a digital «q«nsion stage 12. .s*ere the sampling fi^quency is increased to a significantly high 
integral value, such as a suitable power of 2. followed by an analogue interpolation stage 14 
25 vvhere sample values, at points corresponding to outpm sampling fiequency^ 

Consider a uniformly sampled version x[n] of the band-limited analogue signal x(t). If 
fee samplmg frequency is ^ CTime Period T. ). feen x[n] = xTO- Moreover, if x(t) was band- 
lunited to ^ /2 , then perfect reconstruction can be obtained by applying &e interpolation 
30 function (sampling theorem) 
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x(t) =Ex[k]<j)k(t) where (1) 



n B,^,^ an^l (2) 



Since the summation limit is from -«> to «» it cannot be practically implemented. If non- 
unifoim sampling or finite length is considered (about the point of reconstruction) other types 
1 0 of interpolation functions such as spline, Lagrange interpolator (3) may be used. 



(3) 



The advantage of Lagrange interpolator (3) is tibat it is a polynomial fit constructed in 
15 such a way that each sample is represented by a function which has zero values at all other 
sampling points. 

Evaluating x(t) for all possible values is physically impossible, however there is no 
reason to do so. What is required is to evaluate x(t) at points t = mT, corresponding to re- 
20 sampling with new sampling frequency f = 1/T\ 

Thenn==INT [mP/TJ and ^ = mP - nT, . And therefore, 
x(mT) =x(nT, + At) where 0 < At :s; T 

(4) 
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When converting from a higher to a lower frequency by this method, the output of digital 
interpolation stage 12 is filtered to band-limit the signal content to f72. The filtering process 
removes the effect of aliasing, frequency fold-over due to undersampling. 

5 In converters constructed by figure 1, and op^ting as above mentioned, if the 

interpolator 14 is a Lagrange, spline or linear interpolator, pre-filtering of the extended input 
data generally needs to be effected. Converter 20, shown in figure 2 and constructed m 
accordance with the principles of this invention, depicts a digital e^ansion stage 22 followed 
by sine interpolation 24. By this experience it is possible to avoid said pre-filtering, such as is 

10 necessary when using these Lagrange, spline and linear iateipolators. 

Constraining to the minimum of (t^ , tcF') and limiting the integral of the 
reconstruction formula, equation 5, to [-F/2,F^], sine interpolation may effectively interpolate 
and fiater the expanded data m a single step. In the case where the mput data fireqtiency is less 
1 5 than the output data frequency, the following limit on the integral would apply [-F '/2 JF •/2] . 



F,/2 
•ft/2 



(5) 



Figure 3 illustrates yet another two converters in accordance with the invention. The 
invention has a first stage 32 where the saii5>ling rate of the iiqnit digital data is increased 
20 digitally by an integer &ctor giving the ou^ y[n] at samplii^ rate LF^ The second stage 34 
comprises a simple linear interpolator, vMdi interpolates the denser eaq)anded samples at 
frequency LF, 38 to generate output at required frequency F ' 36. Upsampling reduces tiie 
interpolation error consideiably. Upsampling by a factor of 16 followed by linear interpolation 
leads to SNR of ~ 60 dB. for conversion ratio F '/F, = 4. 

25 

Converts 30 is amplified by nsii^ the same interpoladon fector, 16, for all conversion 
ratios. In effect, the said common interpoladon factor enables ttie same filter coefiBcients to be 
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used for all ratios. A polyphase filter implements the upsampling stage. Figure 5 illustrates this. 

Figure 4 illustrates another converter substantially in accordance with converter 30 and 
in accordance with the invention. The invention has a first stage 42 where the sampling rate of 
5 the input digital data is increased digitally by an integer factor L, giving the output y[n] at 
sampling rate LF^. The second stage 44 comprises a sine interpolator, which interpolates the 
expanded samples at frequency LF^ 48 to generate output at required fi^quency F ' 46. 
Upsampling reduces the interpolation error considerably. 

1 0 Converter 40 is simplified by using the same interpolation factor, 1 6, for all conversion 

ratios, hi effect, the said common interpolation factor enables the same filter coefficients to be 
used for all ratios. A polyphase filter implements the upsampling stage. 

For simple operations converter 30, would be effected in preference to converter 40. 
15 Figure 5 illustrates processes effected in the converters 30 and 40. 

Upsampling, in the embodiments of figures 3, 4 and 5, is generally performed by 
inserting I-l zeros between every two consecutive samples and then filtering the e>q>anded 
result If the converter is constracted in accordance with figure 3, then filtering is performed in 
20 a single step as a part of the interpolation process. 

Insertion of I-l zeros means that Y(z) = X(2^, where y*[n] is the sequence generated by 
inserting I-l zeros in x[n]. In the firequency domain Y'Ce'*) ~ XCe**'), which essentially means 
that the spectrum of x[n] has been compressed I times. Since 'KQ^'^ is periodic in 2n this leads 
25 to creation of extra images in the spectrum. In the case of an embodiment constructed in 
accordance with figure 4, these images need to be removed by a filter with band-limit cOg = tu/I. 

Computational efficiency is obtained in the filter structure above by reducing the large 
FIR filter (h[n]) of length M into a set of smaller filters of length K = M/I, where is I is the 
30 interpolation fector. Since the upsampling process inserts I-l zeros between successive values 
of x[n], only K out of M input values stored in the FIR filter at any time are non-zero. This 
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observation leads to the well-known polyphase filters 

Pfc(n)-h(k + iil) k =0,1,..., I-l- 

n=0, 1, ...,K-1 

5 

The set of I polyphase filers can be ananged as a parallel realisation 62, as shown in 
figure 6, where the output of each filter 64 can be selected by a conunutator 66. 

In the case of linear interpolation, two adjacent polyphase filter outputs are required at 
10 each time. Further reduction in computation is achieved by noting that in the case of linear 
interpolation, not all polyphase filter outputs are used in generating the samples at the output 

In a spedfic example of converter 30, the process of figure 6, rate conversion fix)m 16 
kHz. to 44.1 kHz. The input to the system are samples x[n] at 16 kHz. After upsampling by 16 

15 samples y[n] are produced at frequency 16x16 kHz. = 256 kHz. The output {z[0], z[l], z[2], ..) 
at 441.1 kHz. are interpolations of samples pairs { (y[03,y[l]) , (y[5].y[6]) , (y[ll].y[12]). ..•> 
concsponding to the ratio 256/44.1 = 5.8049. Since only specific points are required others need 
not be computed. In Ihe polyphase approadh, above described, each polyphase filter 62 generates 
an output 64 and the commutator 66 moves to the next polyphase. At the end of a cycle the 

20 commutator returns to the first filter. 

Since only specific polyphase outputs are requked computation can be reduced by 
skipping those polyphase filters whose ou^ut are not required fot that period of time. Unless 
the conversion ratio is an integer no polyphase filter can be absolutely avoided. The above 
25 described example, achieves a computation gain of about four is achieved. 

Internal clock inconsistencies may be a problem in digital fi»quency conversion. 
Consider the example of conversion from 32 kHz. to 44.1 kHz. Real-time systems work on 
limited buffer space and on blocks of data. Suppose tibe constraint on the system is that it always 
30 operates on N output samples. Each time N samples are transmitted at the output the system 
receives an interrupt for DMA (Direct Memory Access) and all the samples collected at input 
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since the last DMA is copied to internal buffer. Similarly N samples must be ready to be 
transferred to the output buffer. 

Now, the input and output clocks are fiee running so there is no guarantee that the ratio 
5 between the time periods of the two clocks will be exactly as computed. As a result it may 
happen that either the number of samples obtained &om input is too few to produce N samples 
at output or they produce more than N samples. 

If Fj is the input sampling jfrequency and F' is the required output sampling frequency, 
10 each time N samples are transmitted at output, [N*Fs/F'] samples should accumulate at the 
input. A small deviation may occur, but on average the above relation must hold. In a case 
where the deviation is appreciable, samples may have to be dropped. This case arises when the 
input rate is hi^er than the output rate. As a result of being dropped samples may have to be 
repeated. 

15 

In accordance with another aspect of the invention, when the input data frequency is 
higher than the output data frequency, more samples are produced at the output, than the buffer 
68 can hold. Overwriting the older samples in the buffer produces a discontinuity. 

20 In the cross fading schraie of figure 7, the ou^ut data frequency 72 is less than the input 

data frequency 71 . To enable a smooth transition, the input data frequency is feded in over time 
as indicated by 73 and flie output data frequency is faded out over time by 74. The result is the 
smooth transition 75. 
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THE CLAIMS DEFINING THE INVENTION ARE AS FOLLOWS: 

1. A method for conversion of input audio frequency data, at an input sample frequency, 
to output audio frequency data, at an output sample frequency, including the steps of: 
5 (a) subjectinglheinputdatatoracpansiontoproduceexpandeddataat^ 

sample frequencj^ and 

(b) interpolating the expanded data to produce the output data, wherein the step of 
interpolating also filters the ou^ut data. 

10 2. A method as claimed in claim 1, where the step of interpolating includes sine 
interpolation substantially in accordance with a reconstruction formula in which an intergral has 
finite limits with respect to frequency to effect filtering during said sine interpolation. 

3, A method as claimed in claim 2, wherdn said smc mterpolation has the constraint that 
15 a cut off frequency, (o^ is equal to pi multiplied by either the input sample frequency, F^, or the 
output sample frequency, F, which ever is lesser, and the reconstruction formula is substantially 
in accordance with: 



20 



4. A method as claimed in claim 2, >?(^erein said sine interpolation has the constraint 
that the cut off frequency, ta^ is equal to pi multiplied by either the mpat sample frequency, 
25 or the ou^ut sample frequency, F, which ever is lesser, and the reconstruction formula is 
substantially in accordance with: 



^' J-F/2 



Ljb=-« 
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5. A method for conversion of digital input audio frequency data, at an input sample 
frequency, to digital output audio firequency data, at an output sample frequency, including: 

(a) upsampling the input data by an mteger factor, so as to increase the sampling 
5 . rate of the input data to produce expanded data; and 

(b) interpolating the expanded data to produce the output data, 

6. A method as claimed in claim 5, where the interpolating is linear interpolatiorL 

10 7. A method as claimed in claim 6, wherein samples for said linear interpolation 
substantially comprise the denser samples of said expanded data, and at least partially 
exclude the less dense samples of said expanded data. 

8. A method as claimed in claim 7, wherein said upsampling includes polyphase filters 
IS for filtering, producing said e3q}anded data. 

9. A method as claimed in claim 8, wherein said polyphase filters are arranged in 
parallel and said i:^sampling includes a commutator for selecting the ouQ}uts of the filters. 

20 1 0. A method £is claimed in claim 9, wherein the commutator at any one time selects 
only two outputs from the polyphase filters. 

11. A method as claimed in claim 5, where the interpolating is sine interpolation. 

25 12. A method of transition of a signal output at one frequency to signal output at another 
frequency, wherein the signal output at said one fi:equency is faded out over a period, and the 
signal ou^ut at said another firequency is faded in over that period, and both are combined to 
produce the signal ou^ut over said period. 
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13. Apparatus for implementing transition of a signal output thereof from a condition at 
vdiich the signal output is derived from a first input signal, of a first frequency, applied on a 
first input to the apparatus, to a condition at which the signal output is derived from a second 
input signal, of a second frequency, appUed to a second input of apparatus, the apparatus 
5 i^vingmeansforfadingoutthefirstmpulsignaloveraperiodandfedingintbesecond 

input signal over said period, and means for comWning these over said period and applying 
the so combined input signals to said signal output 

14. A digital sequential frequency converter for conversion of input audio frequency 
10 data, at an input sample frequency, to output audio frequency data, at an output s^^^^ 

frequency, including: 

(a) means for subjecting the input data to expansion to produce expanded data at 

an output data sample frequency, and 

(b) means for interpolating the expanded data to produce the output data, by an 
15 interpolation process which also filters the output data. 

15. A digital sequential frequency converter as claimed in claim 14. where the 
interpolating is a sine interpolation substantiaUy in accordance with a reconstruction formula 
in which an integral has finite limits, with respect to frequency, to effect filtering during saxd 

20 sine interpolation. 

16. A digital sequential frequency converter as claimed in claim 15. wherein said sine 
interpolation has the constraint that the cut off frequency. a>^ is equal to pi multiplied by 
either the mput sample frequency. or the output sample frequency, P. which ever is 

25 lesser,andthereconstructionfi)rmulaissubstantiallyinaocordancewi1h: 



F,/2 
Fsfl 



2» x[k]e 
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17. A digital sequential firequmcy converter as claimed in claim 15, wherein said sine 
interpolation has the constraint that the cut off frequency, 0)^, is equal to pi multiplied by 
either the input sample frequency, F^, or the output sample frequency, F', which ever is 
lesser, and the reconstruction formula is substantially in accordance with: 



10 18. A digital sequential frequency converter for conversion of digital input audio 
frequency data, at an input sample frequency, to digital output audio frequency data at an 



19. A digital sequential frequency converter as claimed in claim 1 8, where the 
int^olating is linear interpolation. 

20 20. A digital sequential frequency converter as claimed in claim 1 9, wherein samples for 
said linear interpolation substantially comprise the denser samples of said expanded data, 
and at least partially exclude less dense samples of expanded data. 

21. A digital sequential frequency converter as claimed in any one of claims 1 8 to 20, 
25 wherein said upsamphng includes polyphase filters for filtering, producing said expanded 



5 




15 



oulput sample frequency including: 

(a) means for upsampling the input data by an integer factor, so as to increase 
sampling rate of the input data to produce expanded data; and 

(b) means for interpolating the expanded data to produce output data. 



data. 



22. A digital sequential frequency converter as claimed in claim 2 1 , wherein said 
polyphase filters are in parallel and said upsampling includes a conmiutator for selecting the 
30 ou^uts of the filters. 



BNSCOCIO: <WO 0199277A1.L> 



wo 01/99277 



PCT/SGOO/00093 



-15- 

23 . A digital sequential frequency converter as claimed in claim 22, wherein the 
commutator at any one time selects only two outputs from the polyphase filters, 

24. A digital sequential frequency converter as claimed in claim 1 8, where the 
interpolating is sine interpolation. 
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